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We discuss the properties of invariant measures corresponding to iterated 
function systems (IFSs) with place-dependent probabilities and compute their 
Renyi entropies, generalized dimensions, and multifractal spectra. It is shown 
that with certain dynamical systems one can associate the corresponding IFSs 
in such a way that their generalized entropies are equal. This provides a new 
method of computing entropy for some classical and quantum dynamical sys- 
tems. Numerical techniques are based on integration over the fractal measures. 
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In order to characterize quantitatively properties of a given nonlinear system one often uses the 
notion of the dynamical entropy. It describes the asymptotic changes of the system entropy in time. 
Since analytical computing of this quantity is possible only for a limited number of simple models, it 
is important to develop efficient numerical techniques for this purpose. In this article we propose a 
method of computing the dynamical entropy by averaging the static Shannon entropy. The integration 
is performed over a suitably chosen measure, which in the general case displays fractal properties. 



I. INTRODUCTION 



Chaos in a classical dynamical system can be denned by the positiveness of the Kolmogorov-Sinai (KS) dynamical 
entropy jjj . This quantity, characterizing dynamical properties of a system, and defined via an asymptotic limit (time 
tending to infinity) is in general not easy to obtain analytically. On the other hand numerical computing of dynamical 
entropy from time series requires advanced techniques . Also in the quantum case estimating, so called, coherent 
states (CS) quantum entropy 0-0] is not a simple task. In the present paper we propose a method of computing 
dynamical entropy of a system by finding an appropriate iterated function system with the same entropy. 

An iterated function system (IFS) consists of a certain number k of functions Fi, i = 1, . . . , k, which act randomly 
with given probabilities pi,i = 1, . . . ,k. An IFS may therefore be concerned as a combination of deterministic and 
stochastic dynamics. For sufficiently contracting functions one can prove (under some irreducibility conditions) that 
IFS generates a unique invariant measure (see Sect. II). Generically |||| this measure is localized on a fractal set. As 
it was described, e.g., in the elegant book of Barnsley [[[o| IFSs may be used to produce interesting fractal images, 
or to encode and transmit graphics via computer. For the majority of commonly analyzed and applied IFSs the 
probabilities pi are constant. For example, such IFSs have been used to construct multifractal energy spectra of 
certain quantum systems jy] and to investigate the one-dimensional random-field Ising model JlJ] or second order 
phase transitions G3]. On the other hand, with some classical and quantum dynamical systems one can associate 
in a natural way IFSs with place-dependent probabilities [p~4|— p~8|J^ Jig|l . In the present paper such IFSs will be called 
iterated function systems of the second kind, on the analogy of position-dependent gauge transformations ]2(| . 

We estimate the Kolmogorov-Sinai and Renyi dynamical entropies of certain IFSs of the second kind using various 
numerical methods, which can be also applied in the general case. We use similar procedures to analyze the properties 
of the invariant measures of these IFSs and demonstrate their multifractal character. Eventually, we show that one 
can calculate the entropy of certain dynamical systems constructing IFSs with the same entropy. We give several 
examples that illustrate this new method of computing entropy. 

This paper is organized as follows. In the next section the definitions of IFSs of the first and second kind are 
recalled and their basic properties are considered. In Sect. Ill we discuss briefly several methods of analytical and 
numerical computing of dynamical entropy of an IFS. In Sect. IV we study generalized dimensions of measures which 
are invariant under the action of a one-dimensional IFS. Sect. V presents a detailed analysis of a family of IFSs of 
the second kind and their invariant measures. Certain integrals over these measures are calculated. Moreover we 
compute the Renyi entropies, generalized dimensions and multifractal spectra of these measures. In the subsequent 
section we investigate the connection between one-dimensional dynamical systems and IFSs of the second kind. In 
particular, IFSs associated to asymmetric tent map, logistic map, and "hut map" are analyzed and their entropies 
are calculated. We also show how one can apply this method to compute CS-quantum entropy. Concluding remarks 
are contained in the last section. 

In this paper we present only the results and numerical calculations. For the proofs we refer the reader to a 
forthcoming publication. 



II. ITERATED FUNCTION SYSTEMS AND THEIR INVARIANT MEASURES 



An iterated function system (IFS) is specified by k functions transforming a metric space into itself and k place- 
dependent probabilities which characterize the likelihood of choosing a particular map at each step of the evolution 
of the system. Under certain contractivity and irreducibility conditions one can prove the existence of a unique 
attractive invariant measure for an IFS, as well as ergodic and central limit theorems. Miscellaneous results of this 
type have been established since late thirties (some of them have been proved independently by several authors) - see 
for instance pl|- p3| , p^ 15 2^ [33| and references therein. 

In the present paper we study IFSs T = {Fi,pi : i = 1, . . . , k} that fulfil the following (rather strong) assumptions 
which guarantee veracity of the above mentioned theorems: 

General Assumption: 
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(1) X is a compact metric space; 

(2) Fi : X — > X, i = 1, . . . , k are Lipschitz functions with the Lipschitz constants L, < 1; 

(3) pi : X — ► [0, 1], i = 1, . . ., k are Holder continuous functions fulfilling Yli=iPi( x ) — 1 f° r each 

(4) > for every i£l and i — 1, . . . , k. 

Such IFSs are often called hyperbolic. Unless otherwise stated we assume that all IFSs under consideration are 
hyperbolic. 

Let us recall briefly several basic facts on IFSs. 

The IFS T = {Fi,pi : i = 1, . . . , k} generates the following Markov operator V acting on M(X) (the space of all 
probability measures on X): 

{Vv)(B)=y( Pi{\)dv(X), (2.1) 



where v £ M(X) and B is a measurable subset of X. This operator describes the evolution of probability measures 
under the action of T . The related Markov process can be defined in the following way. As a probability space we 
take the code space SI = {I, . . . , k} 11 ^ and we put P x for the probability measure on SI given by 

P x (h, ...,i n ):= P x ({lu £ n : = ijj = 1, . . . , n}) (2.2) 
:= Pil (x)p i3 (F^ (x)) ■ ■ ■ p in (F in _, (F in _ 2 (. . . (F^ (x))))), 

where x £ X, ij = 1, . . . , k,j = 1, . . . , n; n £ IN. Then the formulae 

K¥) = *l(n)(^(n-l)0 ■ ■ W = X (2.3) 

for x £ X,u> £ il,n £ IN define the requested Markov stochastic process on (f2, {P x } x ex)- 

One can show that for an IFS which fulfils our assumption there exists a unique invariant probability measure fi 
satisfying the equation Vfi — /i (the proof of this claim can be found in [[u]]). This measure is attractive, i.e., V n v 
converges weakly to fi for every v £ M(X) if n — > oo or, in other words, J x u dV n v tends to J x u dfi for every 
continuous u : X — > IR. Thus, in order to obtain the exact value of J x u dfi, it is sufficient to find the limit of the 
sequence J x u dV n v for an arbitrary initial measure v. For instance, taking v equal to a Dirac delta measure 5 X for 
some x £ X we obtain the integral of u over the invariant measure /i as the limit of the sequence 

k 

U n ■= ^ P x(h,---,in)u(Xi 1 ,...,i n )> ( 2 - 4 ) 

where Xi lt .„^ n :— F in (F in _ 1 (. . . (F il (x)))). After Barnsley |T(| (see also Q) we call this method of computing 
integrals over the invariant measure deterministic algorithm. To find the integral numerically we can also employ 
the ergodic theorem for IFSs p|,^2 15[ 2(],|3l],|32| . Any initial point x £ X iterated by the IFS generates a random 
sequence (zo = x, Z\, . . . , z n , .77), where Zi := Zf(ui). Then 

^ n— 1 

u(x)dfi(x) — lim — \ u(zi) (2-5) 

v n— >oo n — ' 

A i=0 

with probability one, i.e., except of a set of measure P x zero. Moreover, if u fulfils the Lipschitz condition we can 
evaluate the rate of convergence in the ergodic theorem utilizing the central limit theorem for IFS Jl^,0 . This leads 
to a probabilistic (Monte-Carlo) numerical method of computing integrals over the invariant measure which was called 
random iterated algorithm by Barnsley flG|j . We have successfully applied both techniques to compute numerically 
various integrals, including those necessary to estimate the dynamical entropy of IFS (see Sect. III). 

Finally, let us look at the evolution of densities under the action of IFS. If m is a finite measure on SI and 

Fi (i = 1, . . . , k) are nonsingular (i.e., m(A) = implies m^Ff 1 (A)j = for each measurable A C X), then the IFS 

F generates a Markov operator on the space of densities (with respect to m) on X, also called the Frobenius-Perron 
operator pa]. It is so if, for instance, X is an interval in M and {Fi : i — I, . . . , k} are diffeomorphisms. It follows 



from (2.1) that the Frobenius-Perron operator M associated with T is given in this case by the formula 

m[ 1 ](x) = J2p 1 (k 1 ^))i(f~ 1 (x)) 



dx 



(2.6) 
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where the sum goes over all 1 < i < k such that x € Fi(X), for 7 a density and x € X. 

If probabilities pi are constant then we will say that an IFS is of the first kind. IFSs with place-dependent 
probabilities will be called IFSs of the second kind (they also appear in the literature under the name of learning 
systems). 



III. ENTROPY OF IFS 



Let /j, be the attractive invariant measure for the IFS T = {Fi,pi : i = 1, . . . , k}. We define the probability measure 
P M on the code space ft — {1, ... ,k} by 

P^h, . . . ,i n ) ■■= Pn(W e O : u)(j) = ij, j = 1, . ..,n}) := / P x (h, ■ ■ ■ , i n )dfj,(x) (3.1) 

J x 

for ij = 1, . . . , k, j = 1, . . . , n; n € IV. 

It is easy to show that this measure is invariant with respect to the shift on f2. 
Now we can define the partial entropies as 



k 



H(n):=- 2^ P^(i 1 ,...,i n )lnP fi (i 1 ,...,i n ), (3.2) 
and the relative entropies by 

G(l) := if(l); G(n) = ff(rc) - H(n - 1), for n > 1. (3.3) 

The dynamical entropy of Kolmogorov and Sinai can be extracted from both sequences, i.e., K x = lim, woo G(n) = 
liniji^oo H (n)/n. The usage of relative entropies is often advantageous, since the convergence of H(n)Jn is slow 
(usually as 1/n), while in many cases the sequence G(n) converges to the KS-entropy exponentially fast [p|]3l|. Note 
that the entropy of a stochastic system like an IFS can be defined in several different ways |J7| . Here we are interested 
in the dynamics induced by an IFS in the fc-symbols code space, which leads to the entropy finite and bounded by 
lnfc. 

The concept of dynamical KS-entropy is based on the notion of the Boltzmann-Shannon entropy function which 
can be multifariously generalized |3§f| . In this paper we discuss the two versions of Renyi entropy defined in Q for 
any real q 5^ 1. The first one, often used in the literature, corresponds to the limit of partial entropies: 



1 1 

K„ := limsup In 

n^oc n 1 — q 



[iM<i>- •>*»)]' ■ (3.4) 



ii .... 7 i n — l 

The other one, based on the notion of Renyi conditional entropy [ ^9[ , is defined via relative entropies: 

K q := limsup—!— In [ P^{i x ,. . . ,i n )(p^(i x , . . . ,i n )/P^{h, . . . ,i n - X )\ . (3.5) 



= 1 



For q = both quantities are equal to the topological entropy Kq — ln/c and the KS-entropy is obtained for both 
quantities in the limit q — * 1. On the other hand, in general, the two versions of Renyi entropies are different (see 
Sects. Ill and V). The computation of entropy K q is more straightforward than K q and an analytical treatment is 
possible in some cases, on the other hand, its relation to the thermodynamical formalism seems to be less clear. 

Both definitions give us some form of the Renyi dynamical entropy for the dynamics generated by an IFS in the 
/c-symbols code space and for the specific partition of this space into k rectangles labeled by the first symbol. Note, 
however, that if one defines the (partition independent) Renyi dynamical entropy taking simply the supremum over 
all finite partitions (as for the KS-entropy), this leads to trivial dependence: K q = 00, q < 1; K q = K x , q > 1 p0[ . 
Consequently, from now on, we shall discuss only the Renyi entropy for the above mentioned ^-elements partition. 

For an IFS of the first kind (with constant probabilities pi) both Renyi entropies are equal and can be written down 
explicitly gHH 

K q = K q = ^—\n{p\+pl + ---+pl) 1 (3.6) 
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for q ^ 1. The KS-entropy is obtained by calculating the limit lim g _»i K q , which gives K\ = K\ = — 2i=i Pi ^Pi- 
Observe that for an IFS of the first kind the value of the entropy d oes n ot depend on the character of functions Fj. 

For an IFS of the second kind one cannot directly apply formula ( |3.6| ), since the probabilities are place-dependent. 
A natural generalization for this case is possible |6|J43j, viz., one has to average the Renyi entropy performing an 
integral over the invariant measure \i 



k 

K q =-^—ln[ f2(Pi( x )) 9 M 



(3.7) 



In the limit q — ► 1, corresponding to KS-entropy, this formula gives 

k 



K\ = - I pi (x) ln\pi (x)]dfj,(a 



(3.8) 



Moreover, one can show that the relative entropies converge to the limiting value K q exponentially 143]. Now, to 
compute the entropy, it suffices to apply one of the two methods of calculating the integral over the invariant measure 
of an IFS presented in Sect. II. 

To estimate entropy K q one may consider the modified IFS T q = {Fi,pi(q) : i = 1, . . . , k} with the probabilities 
pi(q) proportional to pi q , that is, given by the formula 



p i (q)tx) = (p i {x)y/Y;(Pi( x )) q - 



(3.9) 



for x € X, i = 1, . . . , k, and q € IR. 

Note that a similar method was used for one-dimensional IFSs of the first kind in . It is easy to prove that the 
IFS T q satisfies our general assumption, and hence, a unique invariant probability measure [i q exists. Then one can 
derive E3] the following inequality 



(l-q)K q > f lnyWaO)*d/i«(aO 



(3.10) 



which provides a lower bound for entropy K q for q < 1, and an upper bound for q > 1. In examples we analyze (see 
Sects. V.C and VLB) this bound is actually very clo se to the exact value of the entropy K q calculated numerically. 
Furthermore, the integral on the right-hand side of ( 3.10| ) can be relatively easily computed (see Sect. II), whereas 
the convergence in (p.4[) seems to be rather slow, namely as rT 1 . 



IV. DIMENSIONS OF INVARIANT MEASURE FOR IFS 



In this section we assume that a one-dimensional (X C M) IFS T = {Fi,pi : i = 1, . . . , fe} is given, where Fi are 
diffeomorphisms fulfilling the general assumption from Sect. II and the following separation condition: intFi(X) n 
intFj(X) — for i ^ j, i, j — l,...,fc, where intFi(X) denotes the interior of the set Fi(X). Our aim is to calculate 
the generalized dimensions D q of the invariant measure for T . These quantities were introduced and analyzed 
by Grassberger, Hentschel, and Procaccia [ p4p5| ] (for more information see p^-|4^|), and Dq is just the Hausdorff 
dimension of the invariant measure. The correlation dimension D2 for certain IFSs has been recently studied by Chin, 
Hunt and Yorke @. 

Let us consider the following pressure function 



P{q,r){x) =:limsup-ln[ Px(h, ■■.,i n ) g \(F in o ... o FiJ (x)\ T 



(4.1) 



i\ ,. . . — 1 



for q > and t e IR. 

In the sequel we assume that the limit in (4.1) does not depend on x and the generalized dimensions D q are given 
by the formula 



G 



D q =^, (4.2) 
where r(q) = r is the only solution of the equation P(q,r) = 0. For the IFS of the first kind this assumption was 



heuristically verified by Halsey et al. |50(] (see also |47|). Moreover Bohr and Rand |5l|,|52j showed that it holds for 
the IFS generated by expanding maps on the interval ("cookie-cutters"). 

To estimate the generalized dimension D q we can use the technique already introduced in the preceding section. 
Namely, we consider the modified IFS T q ^ T — {Fi,pi(q, t) : i = 1, . . . , k} with the probabilities pi(q 7 r) given by 

p i ( q ,r)(x) = {p i (x)r\F((x)\- T /J2(Pj(x)) q \F^)r ■ (4-3) 

3=1 

for x € X, i = 1, . . . , k, q > 0, and r £ 1R. 

Again it is easy to prove that the IFS T q ^ r fulfils our general assumption, and hence admits a unique invariant 
probability measure /i 9 ' T . Then one can show f43[| that the following inequality holds 

(q-l)D q <(q-l)D q , (4.4) 

where (q — l)D q = f is the solution of the equation 



ln^feOr)) 9 \F[{x)Y T d^{x) - 0, (4.5) 
=i 



for q > 0. 

This provides a lower bound for the generalized dimension D q for q < 1, and an upper bound for q > 1. In all 
the cases we study in Sects. V.C and VLB this bound (which can be relatively easily computed) is actually very 
close to the value of the dimension D q calculated numerically. In order to calculate the generalized dimensions D q 
we use the " box-counting" algorithm, which in this case yields better results than the algorithm of Grassberger and 
Procaccia applied to the time series extracted from the IFS. Note that if |.F/(a:)| = L > for alH = 1, . . . , k, 

then the generalized entropies and dimensions are related by a simple formula K q = —D q In L (a relation between 
both quantities for IFSs of the first kind was examined in |l2]]). 

Scaling properties of the invariant measure could be described with the aid of its multifractal spectrum f(a) = 
inf q {aq + (1 — q)D q } (for more information on multifractal spectrum see P, f4l|[4"7| , |53[ ) . 

V. MULTIFRACTALS GENERATED BY IFS OF THE SECOND KIND 

A. Cantor measures 

Let us consider a family of IFSs {X = [0, 1], k = 2; Ft (x) = x/3, F 2 (x) = (x + 2)/3; px(x) = (1 - a) + (2a - 
l)x, P2(x) = a+(l — 2a)x for x £ X}, where a S [0, 1]. It is easy to see that these IFSs fulfil our general assumption 
for a £ (0,1), which guarantees the existence of a unique invariant measures fj, a and veracity of the other results 
mentioned in Sects. II, III, and IV. For a = 1 one can prove the existence of a unique attractive invariant measure as 
well as the ergodic theorem (but not the central limit theorem) using more refined results which may be found in |l5| ] 
or pfij] . On the other hand, for a = 0, the IFS attracts every measure into a linear combination of two Dirac deltas 
localized at points and 1. Hence there exists a whole family of invariant measures {r<5 + (1 — r)8\ : r £ [0, 1]} in 
this case. 

An IFS of the first kind is obtained for a = 1/2, since the probabilities p\{x) = P2(x) = 1/2 do not depend on x. 
The invariant measure fJ-1/2 is spread uniformly over the Cantor set. The generalized fractal dimension is constant 
D q = Da = In 2/ In 3, which implies a singular multifractal spectrum conce ntra ted at ot\ = In 2/ In 3 with /(ai) = oti. 
The generalized Renyi entropy for this IFS can be directly obtained from (3.6). It gives K q = In 2 for all q £ JR. The 



Cantor measure H1/2 can thus be called both uniform (constant generalized dimension) and balanced (constant Renyi 
entropy) @. 

In the case a = 1 the probabilities are place-dependent (pi{x) = x;p2(x) = 1 — x) and define an IFS of the second 
kind. In order to understand the nature of the measure fj,±, let us consider the iterat ions 7„ = Mi(7„_i) of the 



initially uniform density 70 with respect to the Frobenius-Perron operator Al\ given by (2.G) 
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We simplify the notation by introducing the 'box' functions x — > <d a ,b(x) '■= ®(x — a)Q(b — x), with the Heaviside 
function O given by &(y) = 0, for y < 0, and 0(y) = 1, for y > 0. The uniform density in X can thus be written as 
7o = ©o,i- 



Formula (2.6) allows us to obtain for instance the first two iterations of 70 

7i(x) = 9x® i(x) + 9(1 - 35)62^(0:), (5.1) 

and 

72 (ac) = 3 5 x 2 8 ^i (x) + 3 4 (x - 3x 2 )ea ) |(o;) + 3 4 (3a; - 2)(1 - x)G& t i(x) + 3 5 (1 - x) 2 Gi tl (x) (5.2) 
for x £ [0,1]. 

Similarities and differences between densities approximating the standard Cantor measure ^1/2 and the measure /xi 
are displayed in Fig. 1. Due to constant probabilities, in the first case the measure Hi 12 covers uniformly the Cantor 
set (Fig. la, lb and lc). On the other hand, for the IFS of the second kind, the place-dependent probabilities induce 
a highly non-uniform distribution of the measure \i\ (Fig. Id, le and If). For example, in each connected component 
of the support of the density 7„ it can be expressed as a polynomial in x of n-th degree. Note that, if j n achieves its 
maximum at x„, then J n +i(x n ) = 0. One may expect, therefore, that the measure \x\ is multifractal. 

B. Integration over fractal measures 

Let us now calculate the integrals of certain functions u over the invariant measures and Let us find, 
for example, the mean {ua(x) = x) and the mean square (ub{x) = x 2 ) for these measures. Iterating the uniform 
density 70 by the Frobenius-Perron operator Mj/2 we get the sequences of integrals Ua-h = fx u A(x)ln{x)dx = 1/2 
(independently of n) and /7b„^ 2 ' = f x UB(x)j n {x)dx = 3(1— 3~ 2n ~ 2 )/8. Consequently, two integrals in question read 
f x xdfii/ 2 {x) — 1/2 and f x x 2 d^i/ 2 {x) — 3/8. Computing integrals over the measure \x\ it is advantageous to start 
with an initially singular measure. To demonstrate the convergence rate explicitly we take a one parameter family of 
measures consisting of a combination of two delta peaks localized in both ends of the unit interval: n r = r<5o + (l — t)5\, 
where r £ [0,1]. Iterating this measure with respect to the Markov operator V\ given by fl2.l|) we compute the 



r- 

dependent integrals of both functions U An {1) = | [1 + (— |)^] — r (— |) n -> \ and U Bn (1) = \[l + 2{-l) n ]-r{-\) n -> \. 
Both sequences tend to their limits independently of the parameter r, which contributions into the integral decay 
with n as 3 _n . Since both measures and /ii are symmetric with respect to x — 1/2, the first moments (x) are 
equal, however, already the second moments reveal the difference. 

In a similar way an integral of a function over the Cantor set may be expressed as a limit of the sum of 2™ terms 
(multiplied by the appropriate weights), which probe the function on the ends of the intervals composing the Cantor 
set. In some cases this result can be put into a form of an infinite product. For example the characteristic function 



of the uniform Cantor measure \i\i<z is given by |34| (see also 55 1) 



1 x t 







e Ux dn l/2 {x) = e"/ 2 J] cos(-). (5.3) 



n=l 



Due to fast convergence this form is particularly useful for numerical evaluation. In general, computing the integrals 
over multifractal measures generated by IFSs of the second kind one has to rely on nume rica l methods described in 
Sect. II. For IFSs with small number of functions the deterministic algorithm based on (|2.4|) provides more precise 



results than the random iterated algorithm (2.5). The latter seems to be more efficient for IFSs consisting of many 
functions. 

C. Computing of entropies and dimensions via integration over the fractal measures 

The entropy of the IFSs can be expressed as an integral of the Renyi (or Boltzmann-Shannon) entropy function 



over the invariant measure /x a (see (pjl) and (3.8)) for K q , or estimated by the respective integral over the measure 



Hi (see ( 3.10] )) for K q . Note that, comparing the latter case with the former, the natural logarithm interchanges with 



the integral over X. 

Nu merically computed Renyi entropies K q and K q of the measure /xi are displayed on Fig. 2a (however, formula 
( |3.10 ) is valid only for q > in this case). As expected, both entropies depend substantially on the Renyi parameter q, 
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which means that the invariant measure p\ is not balanced. Note that the inflection point of the curve K q is situated 
not at q — but at some negative q c . Making use of the integrals UatX^-) an d Ub u (^) we obtain analytical results 
K2 — (ln3)/2 and K% = (ln2)/2 directly from ( |3.7| ). Renyi entropies allow one to compute the scaling spectra via 
the Legendre transform: g(a) = ini q {aq + (1 — q)K q } and g(a) = mi q {aq + (1 — q)K q } (see Fig. 2b). The common 
maximum of the scaling spectra gives the topological entropy Kq = In 2. Observe that the spectrum g{a) acquires 
also negative values. This does not contradict the interpretation of the scaling spectrum given by Bohr and Rand 
f3lf , which is applicable for g(a). 

Let us recall that the Hausdorff dimension Dq of p\ is the same as for the standard Cantor measure /ii/2 (or any 
other invariant measure p a for a > 0) and equals Dq = In2/ln3 ~ 0.631. The generalized dimensions are given by 
D q = K q /\n3 and can be fairly approximated by D q (see Sect. IV). In Fig. 2c we compare these quantities with 
those obtained by the "box-counting" algorithm and observe that the difference is very small. As expected, the 
generalized dimension decreases with the Renyi parameter q (for example D\ rj 0.47 and D 2 ss 0.41), which confirms 
the multifractal property of the measure p\ (see also Fig. 2d). For this measure we observed that the numerical 
algorithm, providing reliable results for q>l, definitely ceases to work for negative q. 



VI. DYNAMICAL SYSTEM AND IFS 



Let us consider a dynamical system (quantum or classical) endowed with an invariant measure and a partition of 
the phase space. We shall look for an IFS with the entropy equal to the entropy of the dynamical system with respect 
to the given partition. This IFS represents, in a sense, the backward evolution of the system j|3). Having such an IFS 
we could apply formulae for the entropy of IFS given in Sect. Ill and so we would obtain a new method of computing 
the dynamical entropy of the system. We illustrate this procedure on two examples: the Renyi entropy of certain ID 
dynamical systems and the coherent states (CS) entropy of certain quantum systems. 



A. One-dimensional dynamical systems 



Let us consider a piecewise continuously differentiable map / : [0, 1] — > [0, 1]. We assume that there exist subintervals 
Ai (i = 1, . . . , k) such that [0, 1] = (J* =1 A t , f(Ai) = [0, 1], and |/'| > in the interior of A i7 for each i. Let us suppose 
that / admits an absolutely continuous invariant measure p and let us denote its density by p. The partition {^L}f =1 
is generating in this case, i.e., the dynamical entropy with respect to this partition is equal to the dynamical entropy 
of the system. With the map / we can associate an IFS T = {Fi,pi : i = 1, . . . , k} given by 



Fi(x) = f\- A ]{x) 



(6.1) 



and 



Vi{x) 



p(x) 



Fl{x)\ 



(6.2) 



ii). Note that the 



for x S [0, 1] and i = 1, . . . ,k (this is a particular case of the general construction from |17 
functions (-Fi)* =1 are just continuous branches of the inverse of / (see Fig. 3L 

It is well kn own t hat the measure p is also invariant for the IFS T |I^,[l5|,[l8| . Clearly, the generalized entropies 
(given by ("iA) and (3J3)) are in both cases equal, as the probabilities P^(ii, ■ ■ ■ ,i n ) are the same. In general, the 



most difficult stage in t his c onstr uctio n is to show that the IFS T satisfies the assumptions which guarantee the 
truthfulness of formulae (3/7) and 

In the present paper we analyze three examples: asymmetric "tent" map, "igloo" map (better known as the logistic 
map) and "hut" map given by: 

a) Tent map: /(y) = y/r for < y < r and /(y) = (1 — y)/(l — r) for r < y < 1 (where r G (0, 1) is a parameter) 
with the constant invariant density p = 1. Then -^1(2;) = rx, F^ix) = (r — l)x + 1, pi(x) — r, and P2{x) = 1 — r for 
x e [0, 1] (Fig. 3a,b); 

b) Igloo map: f(y) = 4y(l - y) for y e [0, 1] 
for y e [0, 1]. Then F 1 (x) = (1- 

c) Hut map: f{y) 
ye [0,1]. ThenFi(a;) 
x e [0, 1] (Fig 3e,f). 



In this case the invariant density has the form p(y) = \/{n^J y(l — y)) 
VT^A F 2 {x) = (l + VT^A &nd Pl (x) =p 2 (x) = 1/2 for x e [0,1] (Fig. 3c,d); 

(-1 + ^9 - 16 \y - 1/2Q/2 for y € [0, 1]. The invariant density is given by p(y) = y + 1/2 for 
(x 2 + x)/4, F 2 (x) = 1- ((x 2 +x)/4), Pl (x) = (x 2 +x + 2)/8, and p 2 (x) = (6 - x 2 - x)/8 for 







It is easy to show that all the required assumptions are satisfied here and we can use formulae (3.7) and (3.8) to 
calculate the entropy: 

a) Tent map: K q = K q = (ln(r q + (1 - r) 9 )) /(I - q) for q^l, and K 1 = -(r lnr + (1 - r) ln(l - r)); 

b) Igloo map: K q = K q — In 2 for q £ M; 

c) Hut map: K q = ( m (4-*^^fi))/(l - q ) for q ^ 1 and if x = 1/2 + 2 In 2 - (9/8) In 3. 

Note that, for the hut map, one can hardly obtain such an analytical formula for the alternative version of Renyi 
entropy K q . 

Clearly, Kq = In 2 for each of the three maps. In the cases a) and c) D q — 1 for each q. The dependence of D q on 
q in the case b) is presented, e.g., in ||. 

A similar technique can be applied to other classes of ID maps like, for example, "cookie-cutters" introduced by 
Bohr and Rand in |pl] , |52]| . Let us consider, e. g., a repeller given on the unit interval by f(y) = 3y for y £ [0, 2/3] 
and f{y) = 3y — 2 for y £ [2/3, 1], for which typical (with respect to the Lebesgue measure) trajectories eventually 
leave the interval with probability one. Then the measure "uniformly" localized on the Cantor set is the invariant 
measures for this system. The corresponding IFS given by {Fi(x) = x/3, F2(x) — (x + 2)/3} for x £ [0,1] with 
constant probabilities pi — P2 = 1/2 is just the IFS we discussed in Sect. V. 



B. Quantum systems 

In papers we introduced the notion of coherent states (CS) quantum entropy. This quantity may be used to 
characterize chaos in quantum dynamical systems. Out of entire spectrum of Renyi-like quantum entropies K q ||, a 
special meaning may be attached to K\. Namely, the CS-entropy K\ corresponds to the classical KS-entropy. The 
average of K i over the set of all structureless quantum systems, represented by unitary matrices distributed uniformly 
with respect to the Haar measure on U(N), diverges with the matrix size N as In (AT) fjj. This result provides an 
argument in favor of the ubiquity of chaos in classical mechanics (which corresponds to the limit N — ► oo). 

The method of computing the CS-entropy based on the notion of IFS was proposed in ||[l9| (but see also p6| ). 
Again we have shown that one can associate with a quantum system and a partition of the phase space an IFS with 
the same entropy. 

Here we present an exemplary IFS obtained for the family of spin coherent states, the identity evolution operator, 
the quantum number j — 1/2, and the partition of the phase (which is the two-dimensional sphere in this case) into 
two hemispheres (see [|| for details). For this IFS we have: X = [—1,1], Fi(x) = (-3 + 2x)/(6 — 3a;), F 2 (x) = 
(3 + 2x)/(6 + 3x), pi(x) = 1/2 — x/4, and p 2 (x) — 1/2 + x/4 for x £ X. Large contraction coefficient characteristic 
for this IFS ensures fast convergence of integrals performed over the measures approximating corresponding invariant 
measure. It enables us to evaluate numerically the entropy with an enormous precision. For example, the entropy 
K\ » 0.661 314 33271130, being a quantum counterpart of the classical KS-entropy, is evaluated by the deterministic 
algorithm (|2.4D with 14 significant digits. Such precision could be hardly obtained either with random iterated 
algorithm ( |2 . 5| ) or with standard techniques of time series analysis |57],^8| . The quantities characterizing the invariant 
measure of the IFS: a) Renyi entropies: K q , K q ; b) scaling spectra: g(ct), g(a); c) fractal dimensions: D q , D q ; and d) 
multifractal spectrum /(a) are presented in Fig. 4. The common maximum of the scaling spectra gives the topological 
entropy Kq = In 2, while the same curves intersects the bisectrix at the KS-entropy K\. Fractal dimensions D q are 
computed with the aid of the "box-counting" algorithm and compared with the quantities D q defined in Sect. IV. 



VII. CONCLUDING REMARKS 



We have analyzed properties of IFSs with place-dependent probabilities and showed that their invariant measures 
often posses the multifractal property, i.e., the fractal dimension D q depends substantially on the Renyi parameter q. 

We have described a method of computing the generalized entropy for such IFSs by integrating the entropy function 
over their invariant measures. For numerical evaluation of the entropy one can apply the deterministic algorithm 
(useful for small number of functions) or random iterated algorithm (advantageous for large number of functions in 
IFS). Numerical calculations performed for generalized Cantor measures have shown superiority of both methods with 
respect to the standard method of computing entropy from time series generated by IFSs 57,5a]. The entropy and the 
dimension of some IFSs of the second kind studied here display non-trivial scaling properties. The invariant measure 
for such an IFS may be thus neither uniform nor balanced. 
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It is possible to attach an IFS of the second kind to certain dynamical systems in such a way that the generalized 
entropies of their invariant measures are equal. This idea allow us to propose a new method of computing entropy for 
dynamical systems. In this work we demonstrated its usefulness for some classical (Rcnyi-typc entropy of asymmetric 
tent map, logistic map, and "hut" map) and quantum (CS-measurement entropy for two hemispheres, j = 1/2) 
systems. The method of computing entropy by integration over fractal measure has been recently applied to other 
dynamical systems. The tent map with a gap, related to physical problem of communication with chaos, was studied 
in J59| |, while the fractal structure of an exemplary repelling system has been analyzed in p(i|] . Moreover, we used a 
similar method to compute the dynamical entropy of some systems with stochastic perturbations ]6l| . This technique 
may be extended for a wider class of classical and quantum dynamical systems (or even for Markov chains). Such 
results will be presented in a forthcoming publication |ff3| . 

During last few years we enjoyed fruitful collaboration with late Marcin Pozniak. It is also a pleasure to thank Iwo 
Bialynicki-Birula, Lukasz Turski and Daniel Wojcik for inspiring discussions on integration over the fractal measures 
and for indicating the formula ( |5.3p . One of us (K.Z.) is thankful to Ed Ott for hospitality during his stay at the 
University of Maryland and acknowledges the Fulbright Fellowship. Financial support by the Polish Committee of 
Scientific Research under grant No. P03B 060 013 is gratefully acknowledged. 
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FIG. 1. First three iterations of the uniform density on X — [0, 1] by two IFSs {F 1 (x) = x/3, F 2 (x) = (x + 2)/3} attracting 
to the Cantor set. Figures a), b), and c) are obtained for IFS (a = 1/2) with constant probabilities pi(x) = p2(x) = 1/2, while 
figures d), e), and f) for IFS (a — 1) with place-dependent probabilities pi(x) = x, p 2 (x) = 1 — x. 

FIG. 2. Quantities characterizing the invariant measure for "Cantor" IFS of the second kind (a = 1): a) Renyi entropies K q 
(dashed line), K q (solid line); b) scaling spectra g(a) (dashed line), g(a) (solid line); c) fractal dimensions D q (solid line), D q 
(stars); d) multifractal spectrum f(a). 

FIG. 3. Attaching an IFS to a ID dynamical system: a) the tent map, and b) functions Fi and F 2 of the corresponding IFS; 
c) and d) analogous pictures for the igloo (logistic) map; e) and f) analogous pictures for the hut map. 

FIG. 4. Quantities characterizing the invariant measure of the IFS related to the quantum system: a) Renyi entropies K q 
(dashed line), K q (solid line); b) scaling spectra g(a) (dashed line), g(a) (solid line); c) fractal dimensions D q (solid line), D q 
(stars); d) multifractal spectrum f(a). 
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Fig. 4 




